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Abstract 

In this report we conclude the study initiated in Part I and go on to consider optimal 
feedback control (compensator) design for stability augmentation, following the mathematical 
formulation developed in Pan I. We assume co-located (rate) sensors and (force and moment) 
actuators, and allowing for both sensor and actuator noise, formulate stabilization as a 
stochastic regulator problem. Specializing the general theory developed by the author, a 
complete, “closed form’' solution, (believed to be new with this report) is obtained, taking 
advantage of the fact that the inherent structural damping is light. In particular we are able to 
solve in closed form the associated infinite-dimensional steady-state Riccati equations. The 
SCOLF, model involves associated partial differential equations in a single space variable, but 
the compensator design theory developed is far more general since it is given in the abstract 
wave equation formulation. Our results thus hold for any multibody system so long as the 
basic model is linear. 
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A MATHEMATICAL FORMULATION 
OF THE SCOLE CONTROL PROBLEM, PART II: 
OPTIMAL COMPENSATOR DESIGN 


I. Introduction 

In this report we conclude the study initiated in Part I [ 1 1 and go on to consider optimal 
feedback control (compensator) design for stability augmentation, following the mathema- 
tical formulation developed in Part I. We assume co-Iocated (rate) sensors and (force and 
moment) actuators, and allowing for both sensor and actuator noise, formulate stabilization as a 
stochastic regulator problem. Specializing the general theory developed in [2] we are able to 
obtain a complete, “closed form” solution, (believed to be new with this report) taking advan- 
tage of the fact that the inherent structural damping is light. In particular we are able to solve 
in closed form the associated infinite-dimensional steady-state Riccati equations. The SCOLE 
model involves associated partial differential equations in a single space variable, but the 
compensator design theory developed is far more general since it is given in the abstract wave 
equation formulation. Our results thus hold for any multibody system so long as the basic 
model is linear. 

The organization is as follows. To make the report self-contained, Section 2 reviews 
the flexible-structure dynamic equations for a variation of the SCOLE model illustrating in 
particular the generating of the abstract Hilbert space wave-equation formulation. The formu- 
lation of the stabilization as a stochastic regulator problem in the abstract wave-equation setting 
is given in Section 3. The optimal compensator transfer function is derived in Section 4 where 
the steady-state Riccati equtions are also solved explicitly. The compensator is shown to 
have a simple, easily realized structure as a series of band-pass filters centered at 
the closed-loop modal frequencies. 


2. Review: Structure Equations 

We shall derive our results using the abstract formulation developed in 1 1 1. To illustrate 
the generality of this formulation and incidentally also make this report self contained as much 
as possible, we shall re-derive the abstract wave equation formulation for a close variation of 
SCOLE where the controls are concentrated at the antenna end, and for simplicity, exclude 
proof mass controllers. In addition, since we are concerned only with mast stabilization rather 
than the antenna .-dewing problem we shall neglect the kinematic nonlinearity term as well. 

We have then the basic equations for the mast displacements // 0 (-), « e (0 and « v (-), 
(see HI): 


pAw 0 (f, s) + El 0 u”"(t,s) = 0 
pA i/ B (r, s) + EI,u^"(t. s) = 0 ' 


0 < t ; 0 < v < L 


p/ v « v (t, .v) - G/ V «"(t, s) « 0 


( 2 . 1 ) 


where the dots represent derivatives with respect to time t and the primes with respect 
to the space variable .?, as in |1J. With the (“clamped”) boundary conditions at the shuttle 
end (.V = 0): 


«o('.0) = «h('-0) = M'.O) = 0 
»o('. 0) = H 0 ’(t, 0) = 0 


( 2 . 2 ) 


For the controls at the antenna end (.v = L) we have: 


£/ 0 w 0 "'(r, L) = m a u 0 (t, L ) + F y (t) 

l.) + £,(/) 


(2.3) 


where £ (•) and £,(•) are applied forces. With denoting the applied control moment 


(3x! vector): 


£/ 0 w"(r, £) 
FI, »e'(/, L) 
G/ V //y(t, L) 


+ I m.j h 


MM 0 


(2.4) 


where w 4 is the angular velocity defined by: 



0) 4 


«e ('> L) 

«y ('. L) 


(2.5) 


and 1 4 denote the moment of inertia at the antenna end. 

It is convenient at this point to make the transition to a more precise the problem formu- 
lation by employing the vocabulary of “abstract” (or “function space") analysis. Thus let X 
denote the class of 3x 1 functions 


"o CO 
CO 
CO 


0 < s < L 


such that Uq (*). r/J (), u£(-), «„'"(•); «e(0. "eO. «e'(-), «e"'C); « v (0. «y(0> UyO). all 

belong to £- 2 ( 0 , /.), and satisfy the boundary conditions: 

« o (0) - «o(0) = 0 

«e(0) = «e(0) = 0 

«y(0) = 0 . 

Introduce the inner product ini defined by: for any two elements/, g in X, where 



/o 

CO 


So CO 

/CO - 

/e 

CO 

SCO = 

s« (0 


,/y 

(0 


Sy CO 


[/. si = J (/ 0 COs 0 CO + ./i>C0s«C0 + /y COsyCO ) cIs + /o(£-)y’ 0 (C) 

a 

- f ti (L)g H (l.) -< / v (C)s v (£0 *■ /J(Z.) a- 0 '(/.) + /,;(/-).t:,i(M • (2.6) 

The space completed using this inner product will be designated K. It is not difficult to see 
that 

H = Z. 2 [0, Cl 3 x R 5 . 
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It is convenient to use the notation for any element * in H: 


such that for x in X 


X 


( ) 
W H ( ) 
« V () 

bx 

b 2 

h 

b 4 

/>5 


/», = » 0 (L) 
h = i'y(X) 
*>j = u;{L) 
bi = (o 


(2.7) 


(2.8) 


(L) • 

Note that the boundary (at L where the controls are) is part of the state. Define now the 
linear operator A mapping X into 7f by: 


«0 (•) 



«fl (•) 


£7y «„""(■) 

« v () 


-c/ v »;'() 

MO 

; Ar = 

-El 0 u”\L) 

M Z -) 


-Ewm 

"o(L) 


El 0 w o (O 

(O 


f/ 9 »,'(£.) 

« v (O 




Thus defined, A is self-adjoint and nonnegative definite with domain dense in K: 

[.4a , v] = |.v, Ay\ for x , v in X . 

|Av. a j > 0 for a in 1 . 

Then as in [1| we can combine (2.1), (2.2), (2.3), (2.4) to yield the abstract “wave equation” 
in K: 


M\(t) -i ,4.v(0 -» Bu{t) - 0 


( 2 . 10 ) 
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where u(t) represents the control: 


u{t) - 


F y (t) 

FA') 

MAO 


and B is the linear (finite-dimensional) operator mapping R- into K by: 


(2.11) 


Bu - x ; 


0 

0 

0 

u 


(2.12) 


and M is the "inertia” operator: 


Mx = y ; y 


pA« 0 
pA« e 
P ^ ^ 

m a b 2 

b 3 
/ 4 b 4 

bs 


(2.13) 


Note that M is positive definite and has a bounded 


inverse. It is convenient to use the 


notation: 


P Au o 
P Au 6 
P l y « V 
M b h 


( 2 . 14 ) 


where M b is thus a 5x5 matrix, nonsingular and positive definite. In particular, we have for 
the adjoint B * mapping H into /? 5 : 


B*x = b 


B*Mx = M h b . 


We shall call b the “boundary” vector or boundary “trace.” 

Next we need to take into account the random noise input associated with the control. Let 
N s {-) denote white Gaussian noise with (diagonal) spectral density matrix 
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4 ' 

where / denotes the 5x5 Identity matrix. Then we modify (2.10) to read 

Mx(t) + Ax(t) + Bu(t) + BN s (t) = 0. (2.15) 

We are assuming co-Iocated rate sensors. Let v(t) denote the sensor ouptut. Then we can 
model it as: 

v(r) = B*x(t ) + N 0 (t) (2.16) 

where A'o(-) is white Gaussian with spectral density 

d 0 I ■ 

The equations (2.15), (2.16) together yield the state space or abstract formulation that 
we shall employ from now on. 

In concluding this section, let us note the essential properties of A : (i) the resolvent of A 
is compact so that A has a pure point spectrum (ii) zero is in the resolvent set of A — we 
consider only the ‘flexible body” modes. 

3. Formulation of Stochastic Regulator Problem 

The stabilization (or stability augmentation) problem is to make (if possible) the 
“boundary” displacements and displacement rates zero using controls which must be based on 
avilable sensor data. Here the specific version of the problem we shall consider is to minimize 
the steady state “boundary” rates: minimize in other words the time average 

l T 

Mm IWOII 2 dt (3.1) 

7-4* o 

where 

h{t) = B*x(t) , 
b(t) m B*x(t) . 

We shall also impose a soft constraint on the controls and thus consider the problem of 
minimizing 



- 1 - 


1 im 

T 


r T T 

f J ll^(0ll 2 dt + ff IKOII 2 dt 
o 1 o - 


(3.2) 


for fixed \ > 0. This can be seen to be a special case of the stochastic regulator problem in a 
Hilbert space setting treated in |2|. For this purpose we first rewrite (2.15) as a first order (in 
time) equation. Define the state 


Y(t) 


x(t) 

i(0 


(3.3) 


Y( ) is thus an element of the cross product space ftxft. On this space we define, following 
(1), a new norm — the “energy norm — defined on the cross product space: 

T<(vA) x X , denoting domain of ^A , (3.4) 

IIFlj 2 = (svty,, VJy,) + [Afy 2 ,y 2 ] (3.5) 

where 

Y = 

>2 

and 


y, e domain of \s4 . 

Here \A is the unique positive (self-adjoint and nonnegative definite) square root of A, and as 
is well known the domain of \A includes that of A. The cross product space is actually 
complete under this norm, because zero belongs to the reolvent set of A and hence also to the 
resolvent set ol \A. We shall denote the space by The requirement that the first 
component should belong to the domain of 4a will cause no problem for us in what follows. 
Note that calculation of v A is possible using the Balakrishnan formula [6] but will involve the 
boundary “trace" and can be tedious. 

We can now proceed to the reformulation of (2.15) as a Cauchy problem in “state space” 

form in ft (see []]). Thus let 

£ » 

0 / 

-M~ X A 0 


A 


(3.6) 
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25u = 


-M ~ 1 Bu 


Then with 


Y(t) = 


x(t) 


(2.15) can be written: 


Y{t) = JY(t) + 25 u(t) + m v (0 


w'here N s (t) is white Gaussian noise with spectral density 


dj ■ 


It is easily verified that 


where 


Also: 


25*K = -B*y 2 


Y = 


0 -/ 


M A 0 


Let 


C = -SB* 


Then we can write (2.16) in the form 


v(0 = CK(r) + N 0 (t) . 


Finally, the criterion (3.2) now becomes 

r 


I i m 
7 


||M(f)ll 2 dt 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


In this formulation, we may now invoke Theorem 6.9.1 in [2|, where we need to verify 
the stabilizability requirement imposed therein. For this purpose we prove first: 



Theorem 3.1 . (A ~ ®) is Controllable. 


Proof. We need the following Lemma. 


Lemma. Let {<(>*} denote the eigenvectors of A corresponding to the eigenvalues {X k }; 

Wk - , (3.15) 

W4>*. ♦*] - 1 • 

Then the corresponding boundary' trace 


B*$ k * 0 for any k . (3.16) 


Proof. Suppose 


B*$ k = 0 

for some k . Then 


B*(A%) = 0 . 


Since 


we have: 


A$ k = hm k 


Let 


Then 


B*(A$ k ) = 0 . 


u 

u 


0 

0 




A 


«,(« = « p(0 = <(0 = = o 

u e (L) = « e '(L) = u’i{L) = « e '"(L) = 0 
« v (0 = «v(0 = o 
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and of course 

«o( 0 ) = w <?(°) = 0 

«„(0) = « e '(0) = 0. 

In addition: 

£A> w o "( s ) = K u t>( s ) > 0<s<L 

M e C- 5 ) “ w 0 C-^) > 0 < s < L 

Cf/y Wy (j) = Wy (.v) , 0 < s < Lt . 

Since the conditions include those of a beam clamped at both ends, it follows that «<>(•)> u e (-) 
must be identically zero. Similarly the boundary conditions 


Uy(L) = 0 


along with 


« v (°) = M L ) = 0 


are sufficient to make w v ( ) identically zero. Hence B*$ k cannot be zero, for any k. 

Getting back now to the Theorem, let I(t), t > 0, denote the semigroup generated by A. 
Then if (A ~ 35) is not controllable, there must be a nonzero element, say Z, in such that 

[Z, l{t)%u} E =0, r > 0 (3.17) 


for every u in 7?\ Because A has a compact resolvent, so does A and in fact (see also [41): 


V* = 



(3.18) 


V* 


-'“A 


f-L-1 

s/2o )^ ) 


where 


(3.19) 


are orthonormalized eigenvectors of A. Let P k denote the projection operator corresponding 
to the space spanned by \|/*, y k . Then 
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•//OF* = 

•W* - . 


Also, the provide an orthonormal basis in H E and 

2<f)2u = 'LUt)P k $u. 
k 

We note that {vy*. \j7*} is an orthonormal basis for P k K E . Since 


it follows that 




0 


P k %u 


V2 


[u, b k ] {y k ~F*) 


where 


b k - £*<)>* 

Hence 

•UQP&u = -= [«, - «"'"*'?*] = [«,**] 

\ — 

Hence 

[Z, = -\2 [w, ^l(fAf<j, A , 2 ]]o^ sin co t r + [M§ k , z 2 ]co* cos (0*r) 

where 



It follows from (3.17) and (3.24) that 

0 = (m, b k |( 2 |, Mfy k 1 = [w, b k ||z 2 , M§ k 1 , for every k 

Since b k ±- 0 for any k, and u is arbitrary, it follows that 


->/2 (sin co* t )<J> t 
-v2(Ol(cos (O k t)<\> 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


(3.24) 


0 = [zj, M§ k ] - [z 2 , A/<>*] 
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or that 


Z - 0 


which is a contradiction. 


Remark. We may note that while 


we do have that 


B* 4>* * 0 


for any k. 




as k -» oo 


In fact, from 

l|S«l£ = f(lfS«,V*] £ | 2 + |[3w,yj £ | 2 ) 

where the left side 


= [ M b l u , u ] 


and the right side 

= £ II^ShII 2 =f\[u,b k ]\ 2 . 
i '■ i 

It follows in particular that 

£ II** II 2 < °° . 

i 

The result (3.25) has been explored numerically by Taylor and Naidu in |3|. 


Theorem 3.2. Let L b {t ), t > 0, denote the semigroup generated by 

(A - m*) . 


Then 


ll^'b(0^ll f -* 0 as r -t oo 


(3.25) 


(3.26) 


(3.27) 


for every Y. In other words the semigroup l b (t) is strongly stable. 
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Proof. Follows from the general result of this kind due to Benchimol (see [2]) as a 
consequence of the fact that A has a compact reolvent, (A ~ SB) is controllable and A is 
dissipative. 


Theorem 3.2 is adequate to satisfy the stabilizability conditions requirement of 
Theorem 6.9.1 in |2], Applying that Theorem, we note that the optimal control that minimizes 
(3.14) is given by 


«o(0 = - 


SB *P C Y(t) 

X 


(3.28) 


where 


Y(t) 


a - p f m* 


SBSB *P, 


no - p / jo , sv(o 


and 


(3.29) 


K(0) can be chosen arbitrarily. 

Here P c is the unique self-adjoint solution of algebraic Riccati equation 

0 - \P c Y,AY\ e + \AY,P C Y ) F + [SB*y,SB*T] - [SB*P C Y, SB*P C Y] , (3.30) 

Y e 1(A) 

and similarly, Pj is the unique self-adjoint solution of 

0 = [P f Y,A*Y] E + [A*Y,P f Y] E + [d s %*Y,$*Y] - \%*P f Y, d~ 0 '%*P f Y] (3.31) 

Y e SC (A*) . 


The corresponding (minimal) value of (3.14) is 

Tr. 


Tr. SBPySB + 


J-LsLLL 


(3.32) 


We are fortunate that we can solve (3.30) and (3.31) exactly. In fact it is readily 


verified that 
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P e - 'Jxi 


and that 


P f = / • 


Hence 


and 


where 


w 0 (0 


$*f(t) 


V(t) = (A-y < m*)Y(‘) - 



Y 

Correspondingly, (3.32) becomes: 


<dJo + 


x 
V5C ■ 


(3.33) 


(3.34) 


(3.35) 


o + ylld s . (3.36) 

We note that the semigroup i- y (f), t > 0, generated by - ySS*) is strongly stable. We have 
thus obtained a complete closed-form solution of the stochastic regulator problem considered 
herein. 


4. Optimal Compensator Transfer Function 

In this section we shall study the problem of determining the compensator transfer 
function corresponding to the optimal control law (3.33). Strictly speaking (3.34) holds only in 
the weak sense (see [2]), since v(t) contains white noise and is not differentiable. Hence we 
must replace (3.34) by 

Y(t) = j\(r-o)[- V^-)sv(o) da + l/t)Y(0) . (4.1) 

We may set f(0) = 0. Then we have: 

4= J S*l y (t " °)® v (°) d ° . 


«o(0 = 


(4.2) 
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defming the compensator input-output relation. Let 

W(t) = $*i, y (r)S, r > 0 (4.3) 

denote the 5x5 matrix weighting function. Then we need to calculate the Laplace transform: 


f e~ sl S*I y (t) 3 dt , Re. s > 0 . (4.4) 

o 

The integral is absolutely convergent for Re. s > 0, since 

!|®*£ Y (r)$|| -> 0 as (4.5) 

because the semigroup Z y (-) is strongly stable. The Laplace transform 

oo 

J dt = (p/ - (jrf - ySS*))" 1 , p > 0 (4.6) 

o 

is usually denoted 

S(p, A - ySS*) 


and called the “resolvent” of the semigroup (see [2]). (A - ySKB*) has a pure point spectrum 
— or eigenvalue spectrum. To calculate the eigenvalues we need to solve: 


(A - ylKB*)'? = XV ; 


<D, 

o 2 


yielding: 


+ A< I), + Xy BB*d^ = 0 
o 2 = JlO, 


(4.7) 


Let {A.*} denote the eigenvalues, k = 1, .... in order of increasing |A.J. Let, correspondingly, 

X 2 k M* P k + AQ> k + XyBB*<P k = 0 (4.8) 


and let us normalize so that 


[A/dV <DJ = 1 . 


(4.9) 


These <J>* should of course be distinguished from those in (3.15)). Let us note that the {d> t } 
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are not orthogonal. We shall call {A.*} the closed-loop mode frequencies and {<!>*} the closed- 
loop mode shapes. The closed-loop mode frequencies and mode shapes are complex valued. 
Let 


c k - Re. (A.*) ; <S k = Im. (X*) . 

Then, taking inner products with respect to 0* in (4.8) we have: 

X 2 k + r4<D A) <t>J + X*yI|5*4>*|| 2 = 0. 
Let 


b k = Tr. <J> k = B*<t> k . 

We of course assume that y is “small” in the sense that 


( 4 . 10 ) 


YUM 2 « o*] . 


( 4 . 11 ) 


Then 

yIIM 2 
“ “ 2 


<3* = <PJ - ~ ^ = V[yM» t> <t> k ) . 


Consistent with the small y assumptions, we may approximate 


\A< p k , <p k ] 

by 

<t>*l = ro; 

: k = B ^k 

by 

b k = 


using the open-loop (undamped) modes. In this way we get a first approximation to the 
eigenvalues X k as 


, , , y"M 2 

K k = (appr.) - — j — + i(0 k . 


(4.12) 


We may also get a corresponding approximation to <P k as a perturbation of <J) fc . We omit the 
details. The main point is that <f> k will be complex with an approximation of the form 


®k “ §k + ‘ C k&k ■ 


( 4 . 13 ) 
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We can now state: 


Lemma 4.1 . (i) b k * 0 for any k 
(ii) ||£j| ->0 as k -too. 


Proof. Let 




_L 

a k 


K% 


(these 'P* should be distinguished from \j/* in (3.18)) where 

= V[/ldv oj + | 2 . 


Then 'Y k is the eigenvector corresponding to X k . Since 

Tll^ll 2 


<** 


we shall use the strong stability property of the semigroup i. y (*) to prove that 

o k * 0 for any k . 


This is immediate from: 


||/ y (0'»'J = e"*' -* 0 as t -too . 

Hence 


a* < 0, ■ 

proving (i). To prove (ii) we may show that the (0 A ) are linearly independent in X by 
examining the function space part of (4.8). Here rather we shall use what we already have. 
Thus since the {T*} corresponding to distinct eigenvalues are linearly independent, we know 
that 


P*'P A || -4 0 


as 
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Now 


113 * 4^11 


ia*i life* li 


«* 


life* n 

/ \Mt> k f <f>7J 


and 

OJ = 05* + a\ - \X k ? . 

Hence 

life* II -» 0. 

Note in particular that the closed-loop modes <!>* approximate the clamped-clamped modes 
as k -» °o, just as the open-loop mode shapes § k do. 

Next let us note that X k is an eigenvalue if X k is, and has eigenvector W k . Since 

3(p; A-yHB*)'V k = — , (4.14) 

M- A >* 

for any Y in , with the representation: 


I (a k ^ k + b k V k ) 


(4.15) 


we have: 


£(|i ; J - y®B*) Y = 


I 

i 


+ iT^J 


(4.16) 


We want to specialize to 


Y = 3m . 


Also we assume that 3 u is in the (closed) subspace spanned by the {T*}, so that (4.15) holds. 
In physical terms, this is equivalent to saying that we need not consider responses in which no 
modes are excited, We can actually prove this by assuming that there is some (negligibly 
small) strictly proportional damping so that the semigroup /(•) is analytic, implying then the 
same for i, y (-) — see [4], 

Since complex-valued functions and vectors are involved, we now fix our notation for 
the inner product: in C 5 : with asterisk denoting conjugate transpose: 
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Since u is real: 


\u, v] = «*v 


3 )u - 


bar denoting conjugate complex. Hence we can write 

1 




I (M. ql'P* + X I U, cJM** 

L i l 


Taking inner products with respect to and 4 'fj respectively we have: 

V 1 ' \ j 

— r 

+ c£ , m ’W = £ If + 

These equations can be submerged into the form 


where are 2x5 matrices: 


and M kj is 2x2 defined by 


X ?-kM kj = a. | Xjbj Xjbj | 


Zj - | Cj cj | 


^*y - t 


We can solve (4.18) in the form: 


m. Vj\ E m. ni £ 
[n*vni £ in.ni £ 


c > c y I - X X k b k X k b k 

J J 1 


l(°i) 


r jk 


for appropriate 2x2 matrices, r jk , arid hence finally we get 


... . i«. 

93(ki, - YjKB*)5Bm = - y X — — K + = 

1 l 1 ^ “ K k 4 - X* J 


( 4 . 17 ) 


( 4 . 18 ) 


(4.20) 
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and hence: 


K b t c t 


a*S(H, A - tKB*)S - Re. I ^ 


■ K 


^"k b k c k 


\ b k c k 


[a*(^ - I*) a*(lt - x t ) J 


where we note that from 


i m\ r = 

h 


^ I [«, cJ'F* + f [w, c k 1 T* 

*- 1 l 


£ ] *o - — ~ 

« J X ([«, X k b k c%u] + fu, X k b k 
E 1 


(4.21) 




we have that: 


if 

2 i 


(K 

CLe 


b k c t + 




a* 




'5x5 


(4.22) 


Hence finally we obtain that the transfer function (Laplace transform) of the compensator has 
the form 


—yj— i 

vT V <*> 2 


y° ^-k b k C k ^-k b k C t 

1 (H- I*) + cn - x*). 


(4.23) 


From (4.23) we can readily infer the structure of the compensator: we have a series of band 
pass filters centered at the closed-loop mode frequencies 65* , the Q-factor being determined by 
the damping |o* j. The filter amplitudes decrease to zero as the mode number increases, by 
virtue of (4.22). In the limiting case as Jo*| -» 0, the compensator consists of line-filters 
at the open-loop mode frequences (co t /2tt). 

As a first approximation to the solution of (4.18), we may take: 


c j = 


a, 


(4.24) 


u u* 
b k b k 


yielding for the compensator transfer function: 

_l Ji i f b Jl 
vx m> 2 , [at - x,) T (fl _ i k )j 

and a still further approximation by taking 


(4.25) 
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yielding 

. Oi-X^n -*,) ■ (4 ' 26) 

Let us note that in (4.26) we have an approximation for the compensator transfer function 
(Laplace transform) which is quite general — independent of the particular flexible system 
configuration. Only the traces of the open-loop mode shapes on the control sensor boundary 
locations are required — as well of course as the mode frequencies. In particular we note that 
numerical calculations of {b k } for the SCOLE flight article have been reported by S. Joshi [5], 
We can readily see from (4.23) that the feedback amplitude is inversely proportional to VX: 
the smaller the X the larger the control effort, as we expect. Also it increases with the actuator 
noise spectral density and decreases as the sensor noise increases, again as we should expect. 
But the point is that the precise dependence has been determined in (4.23). 
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